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ABSTRACT 

The Parker instability is considered to play important roles in the evolution of 
the interstellar medium. Most studies on the development of the instability so far 
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Q^f have been based on an initial equilibrium system with a uniform magnetic field. 

However, the Galactic magnetic field possesses a random component in addition 
to the mean uniform component, with comparable strength of the two compo- 
nents. Parker and Jokipii have recently suggested that the random component 
can suppress the growth of small wavelength perturbations. Here, we extend 
their analysis by including gas pressure which was ignored in their work, and 
study the stabilizing effect of the random component in the interstellar gas with 
finite pressure. Following Parker and Jokipii, the magnetic field is modeled as a 
mean azimuthal component, B(z), plus a random radial component, e(z)B(z), 
where e(z) is a random function of height from the equatorial plane. We show 
that for the observationally suggested values of (e 2 ) 1 / 2 , the tension due to the 
random component becomes important, so that the growth of the instability is 
either significantly reduced or completely suppressed. When the instability still 
works, the radial wavenumber of the most unstable mode is found to be zero. 
That is, the instability is reduced to be effectively two-dimensional. We discuss 
briefly the implications of our finding. 



Subject headings: instabilities — ISM: clouds — ISM: magnetic fields — magne- 
tohydrodynamics: MHD 
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1. Introduction 

The magnetostatic equilibrium of the system of the interstellar gas and magnetic field 
under the vertical gravitational field of the Galaxy has been shown to be unstable (Parker 
1966, 1967). The physical mechanism for the instability relies on the fact that a light fluid 
(represented by the magnetic field) supports a heavy fluid (represented by the gas) and 
the configuration tends to overturn. It has similarities to the Rayleigh- Taylor instability, 
when a true light fluid supports a heavy fluid. In the Rayleigh- Taylor instability, the fastest 
growing mode has an infinite perturbation wavenumber. However, taking into account a 
"uni-directional" magnetic field along the azimuthal direction in the Galactic disk, Parker 
(1966) showed that the magnetic tension stabilizes large wavenumber perturbations and 
results in a preferred, finite wavenumber. But when the perturbations along the radial 
direction are allowed, those with an infinite radial wavenumber prevail (Parker 1967). As a 
result, the structures formed by the Parker instability are expected to be elongated, and Kim 
et al. (1998) confirmed it through three-dimensional simulations for the nonlinear evolution 
of the Parker instability. 

The Parker instability in the interstellar medium (ISM) has been thought to be a viable 
mechanism in forming giant molecular clouds (GMCs) in the Galaxy (see, e.g., Appenzeller 
1974; Mouschovias et al. 1974; Blitz & Shu 1980). However, the work of Kim et al. (1998) 
raised negative points on that. In addition to the fact that sheet-like structures with the 
smallest scale in the radial direction are formed, they found that the enhancement factor 
of column density is at most ~ 2. The second density issue is eased by noting that the 
interstellar gas can be further susceptible to the thermal instability, as pointed in Parker & 
Jopikii (2000), followed by the gravitational instability. However, the first structural issue, 
which is the direct result of the infinitesimal radial wavenumber, does not easily go away. 

Several ideas on effects that could suppress the maximally unstable nature of the mode 
with an infinite wavenumber have been suggested. One of them is to invoke a "stochastic 
magnetic field" (Parker & Jopikii 2000), which represents the random component of the 
Galactic magnetic field. Using a field composed of the usual mean component and a trans- 
verse component whose strength is weak and random, they showed that in "cold plasma" 
(without gas pressure) the weak, random component exerts a significant stabilizing effect on 
the perturbations with small transverse wavelengths. The physical mechanism is the follow- 
ing. Although weak, the tension of the transverse component that is incurred by the vertical 
gas motions is strong enough to reinstate the gas. They suggested the possibility of preferred 
modes with finite transverse (radial) wavenumbers. Such modes would result in broadened 
structures, which would resemble more the morphology of the GMCs. Their stochastic field 
model is promising in the sense that i) it is consistent with the turbulent picture of the ISM 
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(see, e.g., Minter & Spangler 1996) and ii) it is supported by the observations of magnetic 
field in our Galaxy and spiral galaxies (see, e.g., Beck et al. 1996; Zweibel & Heiles 1997). 
However, their cold plasma approximation needs to be improved. 

The purpose of this paper is to analyze fully the effects of the random component of the 
magnetic field on the Parker instability in a medium with finite gas pressure. We find rather 
surprising results that the random component either reduces the growth of the instability 
significantly or suppresses it completely. And the most unstable mode has a vanishing radial 
wavenumber. The plan of the paper is as follows. Linear stability analysis is carried out by 
analyzing the dispersion relation in §2. Summary and discussion follow in §3. 



We consider the stability of an equilibrium system where gas is supported by its own 
and magnetic pressures against a "uniform" gravity, g, in the negative z (vertical) direction. 
With realistic gravities different growth rates and wavelengths of unstable modes would 
result (see, e.g., Kim et al. 1997), but they make the analysis much more involved. In 
addition, we expect the qualitative features of the stability wouldn't be affected by details 
of gravity. For the magnetic field configuration, the stochastic model suggested by Parker & 
Jopikii (2000) is adopted. It is composed of a mean component, B(z), in the y (azimuthal) 
direction, and a random component, e(z)B(z), in the x (radial) direction. e(z) is a random 
function of z with zero mean. One assumption made on e(z) is that the correlation length 
is small compared to the vertical scale height of the system. So in the equations below, the 
local average is taken by integrating over z for a vertical scale greater than the correlation 
length of e(z) but smaller than the scale height. Then, the dispersion (e 2 ) is taken as a 
constant, which becomes a free parameter of the analysis. With finite gas pressure, p, the 
magnetohydrostatic equilibrium is governed by 



where p is gas density. Two further assumptions are made, which are usual in the analysis of 
the Parker instability: i) an isothermal equation of state, p = a 2 p, where a s is the isothermal 
speed, and ii) a constant ratio of magnetic to gas pressures, a — (1 + (e 2 ))B 2 / '(8irp). Then 
exponential distributions of density, gas pressure, and magnetic pressure are obtained 



2. Linear Stability Analysis 
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dz 



P+( 1 + ( e2 ))^T =~P9, 



(1) 




(2) 



where the e-folding scale height, H, is given by (1 + a)a 2 s /g. 
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The above equilibrium state is disturbed with an infinitesimal perturbation. The per- 
turbed system is assumed to be isothermal too. Since linearized perturbation equations for 
the case without gas pressure were already derived (Parker & Jopikii 2000), the detailed 
derivation is not repeated here. Instead, a reduced form in terms of velocity perturbations, 
(v x ,v y ,v z ), is written down as follows: 
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Here, va is the Alfven speed, B/y/Anp, which is constant over z. Note that the linearized 
perturbation equations for the cold plasma (Eqs. [12] - [14] in Parker & Jopikii (2000)) are 
recovered from the above equations by i) dropping out the terms with a s and ii) noting that 
the scale height of magnetic field (A of their notation) is twice larger than that of gas (H of 
our notation). 



The normal mode solution takes the following form 



•J x , v y , v z ) = (D x , Dy, D z ) exp (J- 4 ik x x 4 ik y y 4 ik z z 4 -^j 



(6) 



where D x , D y , and D z are constants. Taking H and H/ a s as the normalization units of length 
and time, respectively, the dimensionless growth rate, Q = H/(a s r), and the dimensionless 
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wavenumber, (q x ,q y ,q z ) = H(k x , k y , k z ), are defined. Substituting Eq. (6) into Eqs. (3)-(5) 
and imposing the condition of a non-trivial solution, we get the dispersion relation 

Q 6 + C 4 fi 4 + C 2 Q 2 + C = 0, (7) 

where the coefficients C 4 , C 2 and C are given by 

C 4 = 2aq 2 + (2a + l)(q 2 x + q 2 y + q 2 z + 1/4) + 2a(e 2 )(2g 2 + q 2 y + q 2 + iq z /2), (8) 

C 2 = a(a + l)[ql + 4q 2 (q 2 x + q 2 y + q 2 z )] 

+ a(e 2 ) {aq 2 x + (2a + l)q 2 + 4g 2 [(2a + l)(q 2 + q 2 z ) + (4a + l)q 2 y ] + Aaq 2 y [2(q 2 y + q 2 z ) + iq z /2] } 
+ ^ 2 (e 2 ) 2 ql(q 2 x + q 2 y + q 2 z +tq z /2), (9) 

C = 2a 2 q 2 y [2q 2 y (q 2 x + q 2 y + q 2 z + l/A)-(a + l)(q 2 x + q 2 y )] 

+ 2a 2 (e 2 ) { [(a + 1)<£ + 3(a + l)q 2 x q 2 y + (2a + + 4g 2 (g 2 + + g 2 2 ) [g 2 + a(g 2 + g 2 )] } 
" « 2 (6 2 ) 2 g 2 {ql + 2q 2 y - A(q 2 x + q 2 + q 2 z ) [q 2 x + 2a(q 2 x + q 2 )} } . (10) 

Eq. (7) is a cubic equation of Q 2 with complex coefficients. For the case with vanishing 
vertical wavenumber (q z = 0), all the C coefficients become real, and the dispersion relation 
can be easily solved. For small vertical wavenumbers, the imaginary terms, i(e 2 )q z and 
i(e 2 ) 2 q z in Ct and C 2 , can be still ignored. This trick doesn't affect the marginal condition of 
the stability (Q — 0), since C doesn't contain any imaginary term. Here, we remind readers 
of the definition of a. It is reserved in this paper for the ratio of magnetic to gas pressures, 
whereas it was used for the dispersion of e in Parker & Jopikii (2000). For the dispersion 
(e 2 ) is used in this paper. 

Two limiting cases can be considered, which enable us to check the validity of the above 
relation. The formula with (e 2 ) = reduces to the dispersion relation for the original Parker 
instability (see, e.g., Parker 1967; Shu 1974). The C's without the terms containing (e 2 ) 
match exactly with the coefficients of Eq. (53) in (Shu 1974), after imposing the isothermal 
condition 7=1. The other limiting formula is for the cold plasma with p = (Parker & 
Jopikii 2000). As shown above, our linearized perturbation equations recover those for the 
cold plasma. 

The full stability property can be analyzed by solving the above dispersion relation 
numerically. Fig. 1 shows the stability diagram for a — 1. Equi-f2 2 contours with positive 
values corresponding to unstable modes are plotted on the (g 2 , g 2 ) plane for a few different 
values of (e 2 ) 1 / 2 . q z = has been set. Finite g 2 's reduce the growth rate (see, e.g., Parker 
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1966). Three interesting points can be made: with increasing strength of the random compo- 
nent, i) the domain of the instability in the (g 2 , q 2 ) plane shrinks, ii) the maximum growth 
rate decreases, and iii) the g X)max , which gives the maximum growth rate, decreases and 
reduces to zero eventually. Note that without the random component, (e 2 ) 1 / 2 = 0, the most 
unstable mode has the growth rate Q 2 = 0.172 and the radial wavenumber q XtIDax — > oo. 

The above points can be seen more clearly in Fig. 2, which shows the growth rate and 
two horizontal wavenumbers of the most unstable modes as a function of (e 2 ) 1 / 2 for three 
different values of a. Again, q z = has been set. Note that the scale height H changes with 
a. Hence, both the growth rate and wavenumber in real units scale as 1/(1 + a). Even after 
this factor is taken into account, the maximum growth rate increases with a, due to enhanced 
magnetic buoyancy. Two additional points can be made: i) the critical value {e 2 ) l J 2 , above 
which the Parker instability disappears completely, is independent of a, and it is computed 
as l/\/2 = 0.707 from the dispersion relation, and ii) the value of (e 2 ) 1 / 2 , above which the 
radial wavenumber of the most unstable mode vanishes, decreases with increasing a. These 
are the consequences of different roles of uniform and random magnetic fields. 

On the issue of the formation of GMCs, the most interesting range of the values would be 
(e 2 ) 1//2 ~ 0.1, which would result in q x ~ q y for the most unstable modes (although it would 
depend on a, see Fig. 2). Then, the structures formed as the result of the instability would 
be round, so mimicking GMCs. However, observations suggest a larger random component. 
It is generally quoted that in the Galactic plane 0.5 < (e 2 ) 1 / 2 < 1, with the strength of the 
total magnetic field B ~ 3 — 4//G (see, e.g., Beck et al. 1996; Zweibel & Heiles 1997). 
But others such as Minter & Spangler (1996) suggest somewhat smaller values such as 
(e 2 ) 1 / 2 ~ 1/4 — 1/3. If 0.5 < (e 2 ) 1 / 2 < 1, the instability disappears completely or the growth 
rate reduces significantly by more than 80%. If (e 2 ) 1 / 2 ~ 1/4 - 1/3, the growth rate reduces 
by 30 — 60%. But in any case the radial wavenumber of the most unstable mode shrinks to 
zero (q x = 0). That is, the instability basically becomes two-dimensional in the plane defined 
by the azimuthal and vertical directions. This result is opposite to that of (e 2 ) 1 / 2 = 0, where 
the dominant mode of the instability has vanishing radial wavelength (q x — > oo). 



3. Discussion and Summary 

The Parker instability is induced by the magnetic buoyancy of uniform component of 
magnetic field, while gas pressure and random component exert stabilizing effects. The role 
of gas pressure is mainly exercising pressure force along the uniform component, and set finite 
wavenumbers for the instability along the uniform component direction. On the other hand, 
random component threads rising and sinking slices across the uniform field. Its tension 
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becomes stronger at larger wavenumbers. So the role of the random field is to suppress the 
growth of perturbations with large wavenumbers perpendicular to the uniform field. 

Through the linear stability analysis which includes both gas pressure and random 
magnetic field, we have found that with the observationally favored values for the strength 
of the random component, 0.5 < (e 2 ) 1 / 2 < 1, the tension of the random component that is 
incurred by the vertical gas motions is strong enough that the growth of the instability is 
either significantly reduced or completely suppressed. For smaller values, (e 2 ) 1 ' 2 ~ 1/4 — 1/3, 
which are suggested by others, the Parker instability is still operating but with reduced 
growth rate and vanishing radial wavenumber. With (e 2 ) 1 / 2 ~ 1/4 — 1/3, by taking H = 160 
pc and a s = 6.4 km/s (Falgarone & Lequeux 1973), the growth time scale and the azimuthal 
wavelength of the most unstable mode are 70 — 95 Myrs and ~ 2.2 kpc, respectively. They are 
too large for the Parker instability to be a plausible mechanism for the formation of GMCs. 
But it is known that realistic gravity would reduce both (see, e.g., Kim et al. 1997). 
The more serious obstacle in the context of the GMC formation is the fact that the radial 
wavelength of the most unstable mode is infinity. This indicates the structures formed 
would be elongated, in this case, along the radial direction. But it is not clear whether 
such elongated structures would persist in the stage of the nonlinear development of the 
instability. That should be tested by numerical simulations. 

The work was supported in part by KRF through grant KRF-2000-015-DS0046. We 
thank Dr. R. Jokipii for discussions and Dr. T. W. Jones for comments on the manuscript. 
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Fig. 1.— Equi-fi 2 contours on the (g 2 ,g 2 ) plane for (a) (e 2 ) 1 / 2 = 0.05, (b) (e 2 ) 1 / 2 =0.10, (c) 
(e 2 ) 1 / 2 =0.15. In all plots, a = 1 and g z = are used. 
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Fig. 2. — Maximum growth rate, fi^ ax , and its two horizontal wavenumbers, g 2 max and 
<Zy max, as a function of (e 2 ) 1 / 2 for (a) a=0.5, (b) a = 1.0 and (c) a = 1.5. g 2 = for all the 
cases. Note that the growth rate and wavenumber in real units scale as 1/(1 + a). 



